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Abstract 

We analyse here in LO the physical properties of the Green fnnction solntion for the BFKL 
equation. We show that the solution obeys the orthonormality conditions in the physical 
region and fulhlls the completeness requirements. The unintegrated gluon density is shown 
to consists of a set of few poles with parameters which could be determined by comparison 
with the DIS data of high precision. 
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1 Introduction 


The BFKL equation is particularly well suited for description of the behaviour of gluon 
density in the low-x region. Its application in this region is of major importance for the LHC 
and cosmic ray physics. In recent years we have therefore investigated (Tl|^ the solution of 
this equation using the discrete eigenfunction method, hrst proposed in The method is 
closely connected to the Green function approach which is, in our view, the most suitable 
since it does not require any cutoff on the BFKL dynamics. 


The results reported in these papers where as interesting as they were puzzling. The 
proper description of data was only achieved by involving a large number of eigenfunction, 
0(100), contributing in a slowly convergent way. On the other hand, the third and higher 
eigenfunctions were already sensitive to physics at very large scales, much beyond the 
region of data. If true, this property would offer a framework for investigations of Beyond 
the Standard Model (BSM) physics at large energy scales using relatively moderate beam 
energies combined with precision measurements. 

The slow convergence of the procedure given in 0@ suggests, in particular, that the 
discrete eigenfunctions might not by themselves form a complete set. Therefore, in a recent 
paper we rederived the BFKL Green function using the complete BFKL spectrum and 
showed how the imposition of both UV and IR boundary conditions leads to a discrete set 
of poles. The purpose of this paper is to investigate the properties of the full Green function 
solution and to relate it to our previous work. 


The paper is organized as follows: In sectionj^we recapitulate the results of which were 
obtained in the semi-classical approximation. To illustrate the role of approximations used 
in our method we discuss in section |^the case in which the coupling runs without thresholds. 
In this case it is possible to obtain the Green function in analytic form, which reduces the 
solution to an expression in terms of Airy functions in the diffusion approximation in which 
the characteristic function, is simplihed to a quadratic function of oscillation frequency, 
u. In section]^ we discuss a restriction that is imposed on possible paths for the integration 
over the Mellin transform variable, u, arising from thresholds in the running of the coupling, 
and discuss the agreement of numerical results using two substantially different cj-paths. 
In section we discuss the orthonormality and completeness of the BFKL eigenfunctions 
obtained within the approximation that we are using and report that in order to obtain the 
required completeness relation it is necessary to include the continuum of states for which uj 
is negative. In section we discuss the behaviour of the residues of the poles as the gluon 
transverse momentum increases and show that the leading pole is attenuated so that the sub¬ 
leading poles acquire an ever-increasing significance as the transverse momentum increases. 
In section we present the results for the unintegrated gluon density in a model in which a 
very simple ansatz is used for the proton impact factor and for the phase of the oscillations 
of the wavefunctions in the infrared regime. In section we discuss the convergence of 
the sum over poles and show that, in contrast to the results of ref. [^[^, in the formalism 
that we are using here this convergence is quite rapid. Section 7.2 discusses the role of the 
continuum for negative u: on the calculation of the unintegrated gluon density and section 
7.3 shows the relation to the DGLAP calculation. In section 7.4 we discuss the behaviour of 
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the gluon density as a function of x at fixed low virtualities, k^, and show that qualitatively 
it has properties similar to data and display clear contributions from subleading poles. In 
section we show how the presence of new physics albeit at very large energies affects the 
running of the coupling such that both the positions and the residues of the poles are altered 
and this in turn gives rise to a measurable change in the unintegrated gluon density. Section 
[ 9 ] contains our conclusions. 


2 Green Function of the BFKL equation 


We showed that the (Mellin transform of the) Green function for the BFKL equation 
with running coupling can be solved in the semi-classical approximation in terms of Airy 
functions, Ai{z) and Bi{z), where z is a. function of the transverse momentum of the QCD 
pomeron, kx, and the Mellin transform variable, u. 

More precisely the equation for the Green function is 

u}QLo{t,t') - j dt'',t') = 6{t-t'), (2.1) 


where t = ln(fc|,/A^). We note that we have introduced the running of the coupling in such a 
way that the hermiticity of the kernel is preserved so that its eigenfunctions form a complete 
orthonormal set. Eq.(2.1) has a solution, in the semi-classical approximation, of the form 

Quj{tA') = T^N^{t)N^{t')[Ai{z{t))Bi{z{t')9{t-t')+t ^t']. (2.2) 


The argument z{t) of the Airy functions is given by 


{-z{t)f^‘^ = ^ dxu^{x), (2.3) 

where 

xiMt)) = ^ (2.4) 

as{t) 

and x(z/) are the eigenvalues of the kernel /C. The parameter tc is the (ca-dependent) value 
of t for which = 0. Ni^{t) is a normalization factor given by 


N^t) 




(2.5) 


This Green function is analytic in the entire ca-plane with the exception of an essential 
singularity at a; = 0. 


The expression on the RHS of eq.(2.2) has the desired ultraviolet behaviour, namely it 


is exponentially attenuated as f —)■ 00, but the infrared behaviour (where t is small) is not 


properly determined. To obtain the most general solution to (2.1) this Green function should 
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be generalized by adding to Bi the solution of the homogeneous BFKL equation, i.e. the 
transformation 




and eq.(2.2) becomes 




Bi{z{t)) = Bi{z{t)) + cot {(j){u})) Ai{z{t))., 

[Ai{z{t))Bi{z{t')6{t — F) + t -H- 1'~\ . 


( 2 . 6 ) 

(2.7) 


The transformation (2.6) plays the dual role of providing a set of poles at the values of 
uj for which the function = mr and fixing the phase of the oscillatory behaviour of 

the Green function for very small t (or t'), thereby providing the connection between the 
determination of the infrared behaviour of the Green function and the position of the poles, 
originally suggested in (sl. 


It was pointed out in that upon inverting the Mellin transform of the amplitude, 
by integrating along a suitable path in the ca-plane, the saddle-point approximation, used 
to match the BFKL solution with the result of a DGLAP analysis [^[^|^, could be valid 
provided the saddle-point was to the right of all the poles in the ca-plane. On the other hand, 
if this is not the case then the selected contour for the integral over u must surround the poles 
to the right of the saddle-point and will provide signihcant supplementary contributions to 
the unintegrated gluon density which are not matched to the DGLAP result. 


In this paper we report on a numerical analysis of the above-mentioned semi-classical 
solution and discuss the behaviour of the eigenfunctions of the kernel. We confine ourselves 
to the leading order BFKL kernel. The effects of the large components of the NLO BFKL 
kernel will be discussed in a subsequent paper. 


3 Explicit Solution in the Absence of Thresholds 


We begin by considering the simplified case in which the running coupling is given by 

= A ( 3 . 1 ) 

and 00 is a constant. 

Gonsider the eigenfunctions, fuj{t) of the kernel with running coupling 




(3.2) 


The eigenfunctions may be written in integral form as 

f^{t) = J^ [ di/^^(z/)e*^* 


27ruj 


(3.3) 
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with 

In {gM) 

Provided is square-integrable, the contour C may be taken to be along the real axis in 

the i/-plane and we have 


l3oOj 


xW)dy'. 


(3.4) 




2tiu 


dug^{u)e 


iut 


(3.5) 


These eigenfunctions obey the orthonormality relation 

j dtf^{t)f*,{t) = 27 i6{u-u') 


and completeness relation 

J dufjt)f’{t') = 2nd{t-t') 


(3.6) 


(3.7) 


For the leading order BFKL equation with running coupling, g^^ is given by 


gM =g*u2i-^) = {g*M) ^ 


g2i7Ei//(/3o(j) 


T (I +ih') 
F(|-zz/) 




(3.8) 


The integral over u in eq.(3.3) can be approximated by a Gaussian integral around the 
saddle-point, i^ui(t), given by 

xMt)) = (3.9) 


For sufficiently small t the eigenfunctions have an oscillatory behaviour 

1 


fidit) oc 


t 


TT 


x'Mt)) 


cos - -F tu^it) 


(3oOJ 


dv'xW) 


(3.10) 


In terms of the function guj{i^), the Green function may be written (see [^) as 




Airuj 


dp / dp'[ie{p + p') + c{u)]gu:{p)gu 2 {y'), (3.11) 


' —CXO J —OO 


where c{uj) is an arbitrary function of u and this second term reflects the fact that one can 
add to a Green function any solution to the homogeneous part of the equation for the Green 
function. If we write 

c{u) = cot {4>{uj )), (3.12) 

then in the semi-classical limit, we can interpret (p{oj) as being the phase of the oscillations 
of the BFKL eigenfunctions at some small value of t. 
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That (3.11) is indeed the solution to the equation for the Green function can be seen by 
applying the (hermitian) operator 


and using 


to get 

j dTO{t,T)g^{T, t') 




du 


(3.13) 



(3.14) 


Integrating over z/ by parts and using gujiy—v) = {g^jiyu)) ^ yields 

Vid 


dTO{t,T)g,j{T,t') = 


dvr 


/ ‘'G 

’ —CO J —CO 


du'25{u + gui{i^)guj{j^') = 


(3.15) 


Thus we see that it is the factor e(z/ + v') inside the integration over v and v' that 
generates the required inhomogeneous term in the Green-function equation. 


For small t {t < t') this Green function has an oscillatory t behaviour 

f 71 1 

t') ~ sin ( (j){u) + J dv'xW 


(3.16) 


The phase of these oscillations hxed by the boundary conditions at some small value, toi 
of t determines the function (j){u) and hence the positions of the poles at a; = where 
0 (a;) = nvr. 

The Green function has a spectral representation in terms of these poles, namely 




„ (o; - Un) 


and we see from the completeness relation that 

, 6(t — t' 

hm g,j{t,t) = - 

u]—>-cx OJ 


(3.17) 


(3.18) 


In the limit t' —)■ —oo (keeping t hxed), the integral over u' may be approximated by its 
contributions in the regions of the two saddle points at 


u' = 


6 













where 


xMt)) = Po^t- 


(3.19) 


In this limit, > \t^\ and so the function e(i/ + z/') is replaced by ±1 at the saddle points 

respectively. 


We define 


s^{t) = / x{x)dx. 

PoOJ Jq 


(3.20) 


Performing the gaussian integrals over u and u' around the two saddle-points we obtain 

1 


j- irt/^ 

lim ^ 


+ cot (0(a;)) 


giSujit') 


v^ 





(3.21) 




+ 


fUt) 


=~\j (^t') ^ D ^ ^ D} 

Exploiting the symmetry under t -H- t', we arrive at the semi-classical approximation for 
the Green funcion in the region 

, 41n2 

t,t -c G = — 

PoUJ 


Guj{t,t') = 


TT 


x/W 


v/|x'(w(t)||x'(w(t')l 

X 


= |0(t-t')sin (s^(t) + 0 


cos ( s^{t) + 0 + cot(0(a;)) sin + ^ ) + t -H-1' \, (3.23) 


This is the approximation to the Green function given in eq.(2.7) when t,t' tc- A 


similar argument, using the semi-classical approximation, can be used to show that the 


expression (2.7) matches the semi-classical approximation when t or t' S> tc However, the 


semi-classical approximation breaks down if t oi t' 'ph tc, since Uajit) becomes very small and 
the curvature at the saddle-points becomes small. Nevertheless, as we have shown in i , in 
this limit the characteristic function may be approximated by a function which is quadratic 
in u such that the homogeneous part of the equation for the Green function reduces to Airy’s 


equation and the particular combination of Airy functions given in eq.(2.7) generates the 


appropriate inhomogeneous term, so that the expression (2.7) is a good approximation to 
the Green function in all regions of t and t'. 

Henceforth we take the running coupling to be given by 


0!g(t') tt^(to) 


/3o{t')dt' 


(3.24) 


' *0 
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where /3o{t) has steps at the heavy particle thresholds. With this more realistic function for 
the running coupling, it is no longer possible to solve the Green function analytically, even 
in integral form, and we conhne ourselves to a numerical analysis within the semi-classical 
approximation for which the integral over the variable u has a saddle-point at 


4 Numerical Solution Using Different Paths on the uu- 
Plane 


The unintegrated gluon density g{x, t) is given by the inverse Mellin transform of the Green 
function by 


eb2 r 




(4.1) 


where ‘hp(t) is the proton impact factor. We take the Green function to be given by eq.( |2.7 ) 
and the running coupling given by eq.(3.24). The function (j){u) is given by 


TT 


0(u;) + ^ + / = r]NP 


(4.2) 


'to 


with tc given by the relation 


x(0) = 


U 


idc) 


This means that at some small value, to, of t the phase of the eigenfunction with eigenvalue 
uj is given by a non-perturbative phase, rjpfp. 


For a numerical evaluation of the integral over uj, it would be most efficient to identify 
the saddle-point of the integrand and select a path for uj which passes through that saddle 
point and follow the path of steepest descent. In the case where the saddle point lies to the 
left of any of the identified poles of the Green-function the integral must be supplemented 
by the integral over a contour surrounding all poles to the right of the saddle-point. 


Unfortunately, there are restrictions on the permitted paths in the case where the running 
coupling encounters thresholds. In order to consider complex values of u, we require complex 
values of as{t) and hence complex values of t. In the presence of fermion masses m*, the 
running coupling to leading order is given in terms of its measured value at f f by 



where the function F, given by 


F{x) = y ^ In (1 x) + ^/x^ , 

^ Vnp can be a function of lu but in this paper we take it to be constant, although for a realistic fit to 
data we would expect it to possess some w dependence. Note that does not depend on the choice of 
the infrared scale, to, but the infrared phase rjpip does, 
full discussion of this is found in 

















is multi-valued in the complex plane. The running coupling is therefore only uniquely speci- 
hed in terms of t corresponding to kx covering the complex plane once so that the imaginary 
part of t is restricted to 

— TT < < 71, 


which restricts the imaginary part of dg and consequently the imaginary part of u. In 
particular, the calculation of the argument of the Airy functions requires the identihcation of 
ds{tc), so that the restriction on the range of the imaginary part of dg leads to a corresponding 
restriction on the imaginary part of cn. If the real part of uj is small then the real part of 
a{tc) is small, i.e. the real part of tc is large. The restriction on the allowed range of the 
imaginary part of therefore implies that the imaginary part of uj must also be small - i.e. 
we need to select paths which are very close to the real axis in this region. Furthermore, the 
restriction on the imaginary part pushes the possible paths towards the essential singularity 
at very small uj so that we could perform the integration only to some small value of uJmin 
of around 0.05. 


We have selected several paths whose imaginary part differ substantially for large uj. Two 
of them are shown as an example in Figj^ We have performed the integral of eq.(4.1) along 
all of these paths and hnd negligible difference over a large range of t for f = 1 to f = 17 
(corresponding to transverse momentum kx ^ 2 TeV). 

In principle, this result should be expected, since there are no singularities of the Green 
function off the real axis and therefore one can deform the contour into any contour that 
surrounds the real axis and crosses the real axis to the right of all poles. However, our 
contours are not fully closed although the results, at lower t value^ were independent of the 
value. This means that the missing piece of the paths gave a negligible contribution, 
in this t region. The independence of the results from Umin was a hrst sign that the full 
Green function of eq.(2.1) and eq.(2.6l) behaves differently from the slowly converging sum 


of eigenfunctions constructed in ref. |[^[^. In Section 7.1 we will explain in details why this 
happens. 

In addition, the good agreement of the integrals over the different paths shows the nu¬ 
merical consistency of our Mellin transform integration. 


5 Properties of the Eigenfunctions 


The Green function of the BFKL equation 

Giuit, t') = 7iN^(t)Ni^(t') [Ai{z(t))Bi{z(t')9(t - F) -F t -H- F] . (5.1) 

has poles at a; = where 0(ci;„) = rni. If the infrared phase r^^rp is set to a constant 
(ca-independent) value of —vr/4, the hrst few eigenvalues are shown in Table 1. 

^Throughout this paper we have taken the QCD scale A to be 350 MeV 
^e.g. in the t region corresponding to kx < 100 GeV at x = 10“^ 
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omega path 



Figure 1: Two different paths in the cn-plane which were used to invert the Mellin transform 
of the Green function. 


n 

1 

2 

3 

4 

5 

6 

7 

u 

0.389 

0.207 

0.145 

0.113 

0.093 

0.078 

0.0682 

tc 

8.52 

18.7 

29.2 

38.6 

48.0 

57.2 

66.2 


Table 1: The hrst 7 eigenvalues with the corresponding tc values for rj^p = —vr/d 


Except for the hrst three of these, the eigenvalues are well approximated by 

1 

- 1.9(n+f)’ 

consistent with the estimate found in and [^. 

These correspond to a discrete set of eigenvalues of the BFKL kernel with running cou¬ 
pling whose normalized eigenfunctions are given by 

^ (5.2) 

V (p [^n) 

These functions look superhcially like the eigenfunctions of ref. [^[^; the Airy functions 
were dehned exactly in the standard way and the normalization coefficient was previously 
determined numerically (from the requirement that every eigenfunction should be normalized 
to unity). Here the normalization coefficient, is known analytically and is given by 
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Normalized eigenfunctions 



Figure 2: The first three discrete eigenfunctions, /„ plotted against t, showing that these 
eigenfunctions oscillate through (n — 1) nodes and are then exponentially attenuated. 


eq.(2.5); the first factor in eq.(5.2) arises from the conversion between the continuum and 


discrete eigenfunctions. The continuum ones are normalized to a ^-function in u whereas 
the discrete ones are normalized to a Kronecker 5-function in the eigenfunction number, n, 

i.e. 

j = Snn' (5.3) 

This conversion factor can be obtained from the relation 

nUJrL-\-^^rL 

/ 4>' {u:)du: = tt 

JuJn 

where Aun is the separation of the and {n+ 1)*^ eigenvalues. For large n the eigenvalues 
are closely packed and the separation to the eigenvalues may be approximated by 

A _ vr 


The eigenfunctions start by oscillating for small t with n turning points for t < tc and 
are evanescent for t > tc- The first three such eigenfunctions are plotted in Figj^ 
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Figure 3: The normalization factors (2.5) for the hrst three discrete eigenfunctions. 


In Table 2 we show the numerical results for the orthonormality relation (5.3) evaluated 


for the hrst 7 eigenfunctions of eq.(5.2) 


We see that the use of the semi-classical approximation has only had a very small effect 


on the orthonormality relation (5.3). Furthermore, in order to preserve the validity of the 


perturbative expansion, we cannot integrate over all values of t, but only for t > to. However, 
the eigenfunctions are small at small values of t. To see this we note that the additional 
semi-classical factor, N^{t), given by eq.(2.5) is constant for t ^ tc, where the Airy functions 
alone is a good approximation to the eigenfunction, but for sufficiently small t, the variable 
2 ; and x' both become very insensitive to t so that then has a t dependence 

~ 1/ a/cTs (t) - i.e. it decreases in the infrared region as ds{t) grows. The t dependence of 
these normalization factors for the hrst three eigenfunctions is shown in Figj^ 


From Table 2 we hnd that the orthonormality condition is obeyed to a very high accu¬ 
racy. We consider this as a strong indication that we are using a consistent set of physical 
approximations; in particular the semi-classical approximation and the running of as are 
preserving the hermiticity of the Hamiltonian. 
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n\n' 

1 

2 

3 

4 

5 

6 

7 

1 

1.055 

-0.0301 

0.020 

-0.0158 

0.0124 

-0.0107 

0.0093 

2 

-0.030 

1.011 

-0.008 

0.007 

-0.005 

0.005 

-0.004 

3 

0.020 

-0.009 

1.005 

-0.005 

0.004 

-0.003 

0.003 

4 

0.016 

-0.007 

-0.005 

1.003 

-0.002 

0.002 

-0.002 

5 

0.012 

-0.005 

0.004 

-0.002 

1.001 

-0.003 

0.002 

6 

-0.011 

0.005 

-0.003 

0.002 

-0.003 

1.000 

-0.004 

7 

0.009 

-0.004 

0.003 

-0.002 

0.002 

-0.004 

1.000 


Table 2: The ortho normality relation for the hrst 7 eigenfunctions for r]NP = —7r/4 


For the completeness condition we consider the sum 

N 

E(5-4) 

71=1 

For k'rp = 100 GeV we plot this sum in Figj^for the N = 5,10 and 20 eigenfunctions. 

We see that the sum has convergec^ after 10 - 20 eigenfunctions to a distribution on kx 
which is peaked at kx = k'rp. However, we note that the distribution is very broad. This tells 
us that the discrete eigenfunctions do not form a complete set by themselves. Rather the 
completeness requires that the sum over the discrete eigenfunctions must be supplemented 
by the integral over the continuum of states for which uj takes negative values. 

For negative uj there is no critical transverse momentum, tc, beyond which the eigenfunc¬ 
tions diminish, but have oscillatory behaviour for all t. We can write the negative omega 
eigenfunctions as 


f-nit) 


/2 1 


+ r]Np 


'to 


(5.5) 


Since the ultraviolet boundary condition does not impose a specihc ultraviolet phase, the 
spectrum is continuous]^ Although there is no ultraviolet phase-hxing condition, there can 
be an infrared boundary condition which determines the phase of the oscillations at small t. 
For small positive u the eigenfunctions are very closely spaced and become indistinguishable 
from a continuum. For small negative u, the non-perturbative phase should match its value 
for small positive uj in order to ensure a smooth function as u changes sign. Note that for 
large negative uj it may be the case that the infrared phase is not dehned. An example of a 
mechanism in which this happens is where the infrared behaviour of QCD is simulated by 
an effective gluon mass 11 . Here it is found that at some negative uj = —uJi, there is a 
phase transition below which the infrared phase is not determined. Other possible sources 
of such phase transitions could arise from the restoration of conformal invariance at some 
high-energy scale. 


^the fast convergence of this sum will be discussed below 

®This is analogous to the fact that particle only form bound-states for negative energy. Here the analogue 
of energy is —uj. 
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Figure 4: The completeness sum shown as a function of kx and evaluated at k'j, = 100 
GeV, 

eigenfunctions (black). 


(5.4), for the hrst 5 eigenfunctions (red), hrst 10 eigenfunctions (blue) and hrst 20 
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For sufficiently large (negative) values of cu, is given by 




^ exp 



(5.6) 


and 

X' Mt)) ~ 2exp (^E - > (5-7) 

from which we can see that for large t these negative u eigenfunctions have oscillations whose 
amplitude and frequency increase rapidly as t increases. An example for which a; = — 1 is 
shown in Fig. 



Figure 5: The eigenfunction with u = —1, of eq.(5.5) 


With the inclusion of the continuum states, the completeness relation becomes 


lim 




(5.8) 


n=l 


In Fig. [^we plot this quantity for Umin = —1 and Umin = —2 for t' = 10. We see that as 
Umin becomes more negative the LHS of eq.(5.8) becomes more sharply peaked - tending to 
the required (5-function in the asymptotic limit Umin —t —oo. 


6 Transverse Momentum Dependence of the Residues 
of the Poles 

As t increases from t = to, (where the infrared phase is set), the eigenfunction /„ oscillates 
through (n — 1) nodes before the value of tc{n) for that eigenfunction, with (in leading order) 

, , 4 In 2 , , 

tcin) = - ( 6 . 1 ) 

whereas for values of f > tc{n) it decays exponentially. 
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Figure 6: The completeness relation (5.8) for t' 
included for a; > — 1 (left) and u > —2 (right). 


10 , with the negative oj eigenfunctions 




0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

cio 00 


Figure 7: The Green function as a function of ^e{uj} with u close to the real axis. Left 
panel t = 2,t' = 2, right panel t = 12, F = 2 
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This means that for small t, t < tc(l) the Green function as a function of u has a series 
of poles at a; = a;„ with residues that oscillate with amplitudes that decrease with increasing 
n, reflecting a convergence of the sum over pole contributions. This is shown in the left-hand 
graph of Figj^ where we have taken t' = 2 (and u) taken close to the real axis). 

In the region 

G(l) < t < G(2), 

the residue of the leading pole is attenuated and the leading behaviour is now given by the 
first sub-leading pole, which still has an oscillatory residue. This is shown in the right-hand 
graph of Figj^ for t = 12. We note that for this value of t the residue of the leading pole 
has significantly diminished. Note also that the sign of the residue of the sub-leading pole 
is opposite to the case where t = 2, reflecting the oscillatory nature of the residues of the 
poles for t < tc{n). 

As t increases further, the residues of more and more of the sub-leading poles start to de¬ 
cay and the inverse Mellin transform of the Green-function is dominated by the contribution 
from smaller and smaller values of ujn- 

This is in sharp contrast to the situation in which the coupling is kept fixed and for which 
the inverse Mellin (for large rapidity) transform is always dominated by the region close to 
4asln(2). For the running coupling this is not the case, but the value of t at which the 
particular sub-leading poles dominates depend on the value of the rapidity, Y (or ln(l/x) in 
the case of deep-inelastic scattering, leading to an effective pomeron intercept which depends 
on Y. 


7 Properties of the Unintegrated Gluon Density 


The Green function also has a spectral representation given by eq.(3.17) and so we should 
be able to obtain a good approximation to the unintegrated gluon density by summing over 
the pole contributions from n = 1 to n = n^ax- 


As a qualitative demonstration of the unintegrated gluon density that can be obtained 
from the BFKL Green function, we consider a very simple model in which the non-perturbative 
phase, riNP is set to the constant valueof — The proton impact factor has to be positive 
everywhere and concentrated at the values of < 0{1) GeV. It is usually assumed to be 
of the form 


%(kT) = 


(7.1) 


as discussed in ref. 


The form (7.1) vanishes as for small kx, as required by colour 


transparency and the coefficient b has the interpretation of the average inverse square trans¬ 
verse momentum of partons inside the proton and is therefore of the order of 10 GeV“^. The 
overlap integral of the proton impact factor with the eigenfunction for t > to {to corresponds 


^This phase can, in principle, be to dependent although it must lie within the range — f < yNP < f 
to avoid cross-over between eigenvalues. The freedom to select this phase as a function of w is likely to be 
necessary to get a fit between the BFKL formalism and experimental data. 
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to /ct = 1 GeV) is therefore determined by the value of the amplitude at f = to- This is due 
to the fact that the proton impact factor falls for t > to at a rate which is much faster than 
the oscillation frequency of the eigenfunctions in the region t ~ to- In ref also other 
forms of the proton impact factor were investigated, e.g. with different powers of /c^ in the 
prefactor and/or the exponent. It was found, however, that the £t to data has no sensitivity 
to such alternatives due, again, to a small range of the impact factor in comparison with the 
rate of change of the eigenfunction amplitudes. Therefore, for the purpose of this paper it 
is sufficient to take an impact factor which have support only at f = to- 



Figure 8: The unintegrated gluon density (4.1), plotted as a function of fcr, for x = 10“^ 
and X = 10“^ determined with the path (solid line) and pole (dashed line) evaluation. The 
pole result is increased by 3 percent for visibility. 


With these parameters we plot the unintegrated gluon density by performing the inverse 
Mellin transform given by eq.(4.1) over one of the contours shown in Fig. [^and compare it 
with the result obtained from the summation over the hrst 7 poles. We plot this in Fig. [^for 
two different values of x, namely x = 10“^ and x = 10“^. The pole result was increased by 
a factor of 1.03 for visibility. Without this increase both results would be indistinguishable. 
This perfect agreement was obtained because the path enclosed all the poles used in the 
sum. The integration over the path of Fig. [^was performed down to Umin = 0.065, which is 
between the Umin values of the 7th and 8th pole. We checked this agreement for other values 
of Umin and the corresponding sum of poles and obtained an equally good agreement. 


This perfect agreement is, of course, due to Cauchy’s theorem. This agreement is however 
non-trivial because the path integration is not closed what means that the missing piece of 
the path gives a negligible contribution, which is a hrst sign of a very good convergence 
of both the path and pole computations. In addition we note that both computations are 
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numerically very different so it is therefore a very good check of the computational accuracy. 


Another reason for this very good agreement is the fact that we limited the comparison 
to the region of relatively small virtualities, kx < 100 GeV. At x = 10“^, this region is 
experimentally relatively well accessible at LHC and possible future colliders. At x = 10“^ 
and for smaller x values the experimentally accessible kx region diminishes substantially. 
Notwithstanding this, the subleading poles turn out to have a measurable effect (as we shall 
discuss in subsection 7.4) and this will be significant for the prospect, discussed in section 
[^of the identification of new physics from the shift in positions and residues of sub-leading 
poles. 


In the low kx region, as first indicated in Fig. we have a fast convergence of the gluon 
density as a function of the number of poles used in the summation. In Fig. |^we show the 
gluon density computed with different number of poles. The dashed-dotted line shows the 
leading pole contribution, the dashed line shows the sum of 5 poles, the solid line the sum 
of 10 poles and the dotted line the sum over 15 poles. We see that x = 10“^ the first 10 
poles almost converges in the whole kT region, for x = 10“^ the 10 and 15 pole summation 
are completely indistinguishable. We note also that for transverse momenta below around 
10 GeV for the case x = 10“^ and around 20 GeV for smaller values of x, the leading pole 
provides a reasonable approximation. 


7.1 Convergence of the Sum Over Poles 


We can understand how the pole sum converges by considering the behaviour of the normal¬ 
ized eigenfunctions for large n. In this region the eigenvalues are very small {ujn ~ ^/n) and 
the critical momentum tc is very large (proportional to n). This means that for accessible 


values of t the RHS of eq.(2.4) is very small and eigenfunctions oscillate with approximately 


a fixed frequency, uq, given by 


xM = 0 . 


In this region we may use eq.(3.16) and the fact that the phase of the oscillation is t]np at 
t = to to write (j){oj) as 

1 




TT 

Vnp ~ j ~ ^oio + w— 
4 


-I/Q 


xW)dv' 


so that 


0'(a;)- r 


n 




(7.2) 


(7.3) 


The normalization factor N^{t) given by eq.(2.5) has a numerator factor which 

cancels an identical factor in the denominator of the Airy function Ai for t tc and the 
remaining factor is approximately n independent as we replace by pq in the argument 
of x'. 


The upshot of this is that the two factors of 1 /(f)'{u) (see eq.(5.2)) give rise to a 
convergence of the eigenfunction series at small t (relative to G) like ~ 1/n^. Since G 
increases quickly with the eigenfunction number n this fast convergence always happens in 
the experimentally accessible region of t. 
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-2 


Figure 9: The unintegrated gluon density (4.1), plotted as a function of kx, for a; = 10 
and X = 10“^ determined with the pole method. The dashed-dotted line shows the leading 
pole contribution, the dashed line shows the sum of 5 poles, the solid line the sum of 10 
poles and the dotted line the sum over 15 poles. 
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7.2 The Role of the Continuum for Negative uj 


We have seen in section that the contributions from eigenfunctions with negative eigenval¬ 
ues (i.e. negative oj) are essential in order to provide a complete set of eigenfunctions obeying 
the closure relation ( 5.8[ ). In this section we discuss the effect that the Green function with 
negative uj has on the unintegrated gluon density. 


The contribution to this expression for the unintegrated gluon density from the (positive 
uj) poles vanishes asymptotically with t, reflecting the behaviour described in section 
whereby as t increases the residues of more and more of the poles pass from an oscillatory 
behaviour to an exponentially damped behaviour. 


The contribution, A^<o fo g{x, k) from negative uj is given by 


A^<o5'(a^, k) 


lim 

^min^ — OO 





with given by eq.(5.5) (and t 


2 ln(fc'r/A)). 


(7.4) 



Figure 10: Blue line: 
UJ 0 to UJrain 
taken x = 10“^ and t 


The integrand for the RHS of eq.(7.4). Black Line: The integral from 


We see that the integral has converged well for 
= 12 


= —1. We have 
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At first sight, one would expect this contribution to have a negligible effect on the un¬ 
integrated gluon density owing to the factor of x‘^. However, as can be seen from Fig. [OT 


for values of u just below a; = 0 the integrand of the RHS of eq.(7.4) is still quite large 


and rapidly oscillating, although we can also see that the integral converges to a fairly small 
value by cu = —1. 

We have computed the gluon density including the contribution from negative uj. The 
contribution of positive cn’s was given by the sum of the first 10 poles and the contribution 
between oj = 0 and cuio (where the poles are densely packed) was treated as a continuum. 
The negative omega contribution was evaluated from a; = —2 to a; = 0 using eq. (7.4). Fig. 
11 shows the unintegrated gluon density for x = 10“^, x = 10“^ and x = 10“'^ including 
the contribution from negative u. The dashed line shows the pole contribution (computed 
from the sum of 10 poles), the solid line shows the same pole gluon density with added 
contribution of negative w’s. The negative u contribution was computed assuming that the 
infrared phase for positive and negative u match at a; = 0. 


Fig. 


11 shows that the contribution of negative u may be substantial at x = 10 ^ if 

However, this 


the negative u infrared phase is hxed and substantially different from 0 
contribution diminishes very fast with decreasing x as can be seen in the same figure. 


7.3 Comparison with DGLAP 

We see from Figj^ that for x = 10“^ above kx of around 30 GeV, the unintegrated gluon 
density ceases to rise, whereas for smaller values of x, e.g. x = 10“^, the unintegrated gluon 
density continues to rise up to transverse momenta of 100 GeV. This can be understood 
from the t-dependence of the residues of the poles as seen for example in Figj^ For t = 12 
(corresponding to kx of around 100 GeV), we see that the residue of the sub-leading pole is 
of opposite sign and slightly larger in magnitude that that of the leading pole (whose residue 
is evanescent since t is above tf. for that pole). The contribution of the sub-leading pole is 
suppressed by a factor of 

^ ^ 0.2 

For X = 10“^ the contribution from the sub-leading pole is sufficient at such values of t 
to halt the rise in the unintegrated gluon density, whereas for x = 10“^ it is insufficient. 
Nevertheless, at sufficiently large t the unintegrated gluon density for x = 10“^ will also 
cease to rise and will eventually display oscillations. Above some large t, outside the range 
of t ( kx) range considered in this paper, these oscillation are certainly unphysical because 
they lead to negative gluon density. 

There is no reason a priori why the BFKL amplitude should not display oscillations. The 
inversion of the Mellin transform consists of an integral over u which has greatest support 
at the saddle-point Us- For values of t below tc for this value of u the amplitude displays 
oscillatory behaviour and it is only when t exceeds this critical value that the oscillations 
halt. The unphysical oscillatory behaviour indicates that the solution to the BFKL equation 
is being applied to deep-inelastic scattering outside the kinematic range for which it was 
intended. The application of the gluon scattering amplitude in the Regge regime to the 
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Figure 11: The unintegrated gluon density for x = 10“^, x = 10“^ and x = 10“^ including 
the contribution of negative ca’s. The dashed line shows the pole contribution (computed 
from the sum of 10 poles), the solid line shows the same pole gluon density with added 
contribution of negative cn’s. 
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determination of the unintegrated gluon density identifies the rapidity with ln(l/a;), which 
in LO is only valid provided the rapidity signihcantly exceeds {t — t'). Therefore for x = 10“^ 
and t' conhned to the low t region where the proton impact factor has non-negligible support, 
we would expect the BFKL determination of the unintegrated gluon density to become invalid 
if t is substantially larger than ~ 5 (corresponding to /c^ of order 10 GeV). 


We would expect this limitation on the allowed range of t to be less stringent when NLO 
effects are taken into account. Indeed, as pointed out by Salam 13 , the LO treatment 


ignores the discrepancy between the rapidity variable used in a BFKL analysis (which is 
symmetric in t and t') and the variable x. This introduces a factor of 




whose absorption generates the largest part of the NLO contribution to the characteristic 
function, x{^)- We would therefore expect the BFKL amplitude computed at NLO to be 
less sensitive to the difference {t — f) than a purely LO analysis. 

It is known that at sufficiently large t and sufficient small x, the double logarithm limit 
(DLL) of a BFKL analysis matches that of a DGLAP analysis. In Mellin space, the region 
in which this approximation is valid is given by 

1 3> cj S> as{t). 

In terms of x this translates into the limits 

« 1, ln(i) » 1 

Moreover the match between a DGLAP analysis and a BFKL analysis can only be valid if t 
exceeds tc a.t u ^ ujg, where Ug is the saddle-point for the inversion of the Mellin transform, 
i.e. the region around Ug is where the integrand has its maximum support (assuming that 
this saddle-point lies to the right of all poles). We have seen that we need to have values 
of X smaller than x = 10“^ in order to avoid a signal from the oscillatory part of the 
BFKL eigenfunctions. This means that in the case of a DGLAP analysis, the DLL limit can 
only be reached for very small values of dg{t), i.e. values of t way beyond any reasonable 
experimentally accessible region. For the case of BFKL with running coupling we need to 
go to even smaller values of x before the DLL becomes a reasonable approximation. This is 
because for running coupling the contributions from leading poles are attenuated at large t 
and we need to be at sufficently large rapidity to ensure that these leading poles dominate 
the unintegrated gluon density. The input for a DGLAP fit is the structure function at some 
reference photon invariant mass, Qq. In the case of the discrete BFKL pomeron the input 
would be the proton impact factor and the infrared phase rjjsip. As discussed in section 
this impact factor is expected to be dominated by its value at f = to so that the only free 
parameter associated with this impact factor would be the overall normalization. The only 
other parameter which can be substantially varied is the infrared phase, rj^p ,which should 
be a function of the eigenvalue Un- The infrared phases are generated, as the eigenvalues 
ojn, by the quasi-bound states of gluons inside the proton and therefore should be described 
by a simple parameterization. Our previous experience indicate that two parameters may 
be sufficient in order to generate the rj — uj dependence required to fit data. 
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Because of these very different parameterizations, it is quite likely that there exists an 
overlap kinematic region at low-a; for which data can be equally well described by a con¬ 
ventional DGLAP £t or a fit to the discrete BFKL pomeron. A detailed comparison of the 
discrete BFKL pomeron with data will be performed in the next paper, after the NLO cor¬ 
rections are introduced. In the same time we will discuss the comparison with the DGLAP 
fit. 

7.4 Dependence on x 



X 


Figure 12: The unintegrated gluon density as a function of x determined from the pole con¬ 
tribution only, at various k^. The dashed line shows the leading pole contribution normalized 
to the values of the gluon density at x~‘^, for each k^. 


It is well known from HERA data that the x dependence of F 2 , which is directly connected 
to the gluon density, exhibits a striking Regge type behaviour, i.e. ~ {1/x)^. The parameter 
A is not a constant, it increases logarithmically with (which we set here equal to k'^), see 
e.g. 


14 


Such a behaviour is also a feature of the gluon density obtained from the Green 
function solution of BFKL investigated here, see Fig. The figure shows the unintegrated 


gluon density as a function of x, determined from the pole contribution only, for various 
values of k^. In this log log plot the function (1/x)^ is a straight line, so we see immediately 
that the gluon densities exhibit the same striking linearity as data. The slope A increases 
slightly with k^ which can be seen by comparison of the slope of gluon density (full line) 
with the A;|,-independent slope of the leading pole contribution (dashed line). The leading 
pole contribution is, for each k'^, normalized to the values of the gluon density a.t x = 10“^. 


In Fig. 12 we display the gluon density in the validity region of our solution to BFKL, 
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X 


Figure 13: The unintegrated gluon density as a function of x determined from the pole 
contribution and the negative u contribution, at various The dashed line shows the 
leading pole contribution, normalized in the same way as in Fig. 


At < log(l/a:), which means that kx should be smaller than order of 10 GeV at x = 10“^. 
At larger k^, not shown in this figure, we observe a clear deviation from linearity in the 
region of x between 10“^ and 10“^. A close inspection of the highest k^ line in Fig. 12 shows 
that this effect sets in already at a relatively low fc-r of about 30 GeV. 


In Fig. [T^we show the x dependence of the gluon density which includes the negative 
u contribution (full line) for various fc|,’s. The dashed line shows for comparison the same 
leading pole contribution as in Fig. We see clearly from this plot that the negative u 
contribution substantially affects the linearity in the region of x > 10“^. Since this non- 
negligible contribution from negative u occurs at the relatively low values of k^ considered 
here whereas data strongly indicates linear trajectories in x up to x ~ 10“^, we take this 
as a strong hint that the infrared phase, rj^p for a; < 0 and the precise form of the proton 
impact factor are such that the overlap of the proton impact factor with the negative uj part 
of the Green-function is very small. In view of the fact that the proton impact factor is 
expected only to have signihcant support near to where the infrared phase, //vp, is hxed, 
such suppression of the overlap would occur if the infrared phase were very small in the 
region a; ~ 0. 


In this paper we do not attempt to make any detailed comparison with data because we 
are working here in the LO only and it is well known that the NLO and LO results differ 
substantially in BFKL. We note, however, that the difference between the full gluon density 
and the leading pole contribution seen in Fig. 12 is due to the sub leading poles. A similar 
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change of slope A with , even more prominent then in onr LO calcnlation, is also seen 
in the data 14 . Therefore it is highly possible that the properties of the snbleading poles 


can be well determined from the measnrements of the slope A since in onr solntion there is 
only one free parameter per pole, and the nnmber of contribnting poles is small, dne 

to the fast convergence of their snm. 


8 The Effect of New Physics 


In previons pnblications , we have pointed ont that both the positions of the poles in uj 
and their residnes are sensitive to any new physics which affects the rnnning of the conpling 
and hence the valne of the critical momentnm, tc, provided tc is above the threshold for snch 
new physics. This nniqne effect is dne to the fact that an eigenfnnction of the BFKL kernel. 


(5.2), is a (qnasi) bonnd state of glnons with very different virtnalities, ranging from to tc- 


Even if snch a state is probed at low t valne, mnch below the Beyond Standard Model (BSM) 
threshold, the resnlt is sensitive to the properties of the whole state since the eigenvalne and 
the residne at the probed t is determined by all the glnons of the stat^ 

Formally speaking the sensitivity to BSM thresholds emerge from the fact that the fnnc- 


tion eq. (4.2), which defines the eigenvalnes Un by the reqnirement = nvr, contains 
an integral over the freqnency which ranges from fg to tc- This freqnency, defined by 
eq. (3.19), is strongly sensitive to a snpersymmetry (SUSY) threshold becanse the valne of 
/So changes snbstantially in the SUSY region. 


In Fig. [^we show an example of this in which we plot the Green fnnction as a fnnction 
of a; on a path close to the real axis, for a typical low t valnes of 9, 6.7 and 4.4 corresponding 
to kp of 30, 10 and 3 GeV. We are comparing the Standard Model (SM) with the MSSM 
with a snpersymmetry (SUSY) threshold of 3 TeV. We see that the position of the hrst pole 
is nnaffected becanse the corresponding tc lies mnch below the SUSY threshold of 3 TeV. 
However the positions of the snb-leading poles are shifted to the right becanse their tc’s are 
either close to the 3 TeV threshold, as in the case of the second pole, or mnch above it for the 
rest. It is interesting to observe that the residnes of the non-leading poles oscillates strongly 
in the displayed t region. This may help to disentangle their contribntion since this t region 
is well accessible to high precision measnrements. 


9 Conclusions and Outlook 

We have investigated here the properties of the complete Green fnnction solntion to BFKL 
eqnation in LO approximation, nsing a mixed techniqne of analytic and nnmerical analyses. 
We have shown that this solntion fnlfills the completeness reqnirement and leads to a set 
of eigenfnnctions which are properly normalized and are orthogonal to each other. These 

® The measured variable is closely related to the BFKL variable k^. 

®Since t = ln(fcf./A^), from Table 1 we see that kp which corresponds to tc of the first pole is 25 GeV, of 
the second one is about 4 TeV, of the third one 500 TeV, of the fourth one 54000 TeV, etc. 
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mathematical properties are fulfilled with a high numerical precision, which is not trivial in 
view of the fact that the eigenfunction states are defined only for t values above some small, 
but perturbative to, and not in the whole region — cx) < f < oo, as would be mathematically 
required. To achieve the completeness it is mandatory to take into account the contribution 
from the states of the negative u continuum. 

The unintegrated gluon density is defined, in this paper, as the inverse Mellin transform 
for the Green function of the BFKL equation over a suitable path in the complex uj plane. 
We show that at low t the integration over an uj path is exactly equivalent to the sum of 
poles supplemented by a possible contribution of the negative uj continuum. The sum of 
poles is dominated by a relatively few terms, which is in contrast to the situation found 
in the discrete eigenfunction solution of [T]-[^. The fast convergence of the complete Green 
function solution is due to the presence of the t dependent normalization factor which was 
missing in the pure eigenfunction approach of m- 

We have investigated the region of validity of the unintegrated gluon density obtained 
in this paper and found that it is limited to the region At = t — t' < log(l/x). If At 
is sizably larger the unintegrated gluon density starts to exhibit an oscillatory behaviour, 
which eventually leads to an unphysical, i.e. negative gluon density. Since in DIS t' is 
limited by the proton factor, which confines it to very small values, the t values cannot be 
too large. At x = 10“^ they corresponds to kx values of the order of 10 GeV. At smaller x, 
like X = 10“^, the virtuality kx would be an order or more of magnitude larger, however at 
HERA or LHG the region of accessible virtuality decreases with x substantially, so that the 
region of applicability of the BFKL equation is limited effectively to kx of the order of 10 
GeV. 


At low kx the gluon density is dominated by the leading pole, which leads to a linear 
behaviour in the logarithmic x dependence with a slope which varies with kx- This variation 
is due to the contribution of the subleading poles. We show here that the uj values of the 
subleading poles are sensitive to the Beyond Standard Model (BSM) effects due to the same 
mechanism of threshold sensitivity as investigated in [^|^. Therefore, the deviations from 
the leading pole behaviour are sensitive to the BSM effects and could be measurable due 
to the high precision of the measured slopes of the logarithmic x distribution, especially in 
future projects, see e.g. 15,16 . However, in this paper we did not attempt to perform any 


data analysis because we are working here in the LO order approximation only and it is well 
known that the u: values differ substantially between LO and NLO approximation of BFKL. 
The LO analysis presented here allows the full understanding of the qualitative properties 
of the BFKL solution in the low kx region which we expect to be the same as in the NLO 
analysis. The quantitative results in the NLO approximation will be presented in our next 
paper in which we will also perform a comparison with the high precision DIS data. 
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